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We identify a new parameter that controls the localiza- 
tion length in a driven quantum system. This parameter re- 
sults from an additional quantum degree of freedom. The 
center-of-mass motion of a two-level ion stored in a Paul trap 
and interacting with a standing wave laser field exhibits this 
phenomenon. We also discuss the influence of spontaneous 
emission. 

PACS numbers: 42.50.-p, 05.45. +b, 03.65.Sq 



I. INTRODUCTION 

The phenomenon of localization manifests itself in 
many quantum mechanical systems ranging from the lo- 
calization of light in a random medium pj, to Ander- 
son localization of an electronic wave Q and to atoms in 
time-dependent laser fields (|^] . In all these cases the un- 
derlying classical system is chaotic and shows diffusion as 
a function of time. In contrast, the quantum mechanical 
counterpart has a localized wave function whose width is 
governed by the classical diffusion and Planck's constant 
I n the present paper we show that there exists an 
additional quantum parameter that controls the localiza- 
tion length. In the system of a two-level ion stored in a 
Paul trap and interacting with a standing wave it is 
the detuning between the transition frequency of the ion 
and the laser frequency. 

A recent paper Q has shown that a stored ion mov- 
ing along a far-detuned standing wave shows localiza- 
tion in position and momentum variables. In this work 
we have neglected the internal structure of the ion. We 
have therefore only considered the quantum dynamics of 
a particle moving in a one-dimensional time-dependent 
potential. In the present paper we extend this analysis 
and consider the motion of an ion taking into account its 
internal dynamics. 

Our paper is organized as follows: In Sec. || we first 
summarize the essential ingredients of the problem. In 
particular, we introduce the relevant equations and de- 
scribe the methods of solution for the classical and quan- 
tum mechanical equations of motion. We also define the 
quantities which we calculate such as the position and 
momentum distributions and their moments. 



In Sec. [Ilj we discuss the dependence of the classical 
and quantum mechanical position and momentum distri- 
butions on the detuning. We find characteristic oscilla- 
tions in the widths of these distributions as a function of 
the detuning. These oscillations are absent in the corre- 



sponding classical curves. We explain these structures by 
transforming the Hamiltonian into an interaction picture. 

We emphasize that this is different to the case of atoms 
moving in a phase modulated standing wave. There the 
corresponding oscillations appear in the widths of the 
quantum as well as of the classical momentum distribu- 
tions pl] |. Hence, the detuning is an additional param- 
eter controlling the width and the shape of the quantum 
distributions. 

Moreover, in the quantum systems discussed so far in 
the context of dynamical localization the quantum dif- 
fusion is always slower than the classical one. However, 
in the present model the new control parameter detun- 
ing can create situations in which the quantum diffusion 
temporarily exceeds the classical one. 

Quantum interference effects such as dynamical local- 
ization are extremely sensitive to decoherence arising for 
example from noise or spontaneous emission |l2[ . We are 
therefore forced to investigate in Sec. IV the influence of 
spontaneous emission. We find that the drastic difference 
between quantum and classical behavior is preserved in 
the limit of a far off-resonance situation. We conclude in 
Sec. [\] by summarizing our main results. 



II. FORMULATION OF THE PROBLEM 

We consider the standard Paul trap set-up realized 
experimentally in many labs p3[-p5|: a standing elec- 
tromagnetic wave of frequency lul and wave number k 
aligned along the x-axis couples the internal states \g) 
and |e) of a single two- level ion of mass m to the center- 
of-mass motion. The resulting dynamics of the ion fol- 
lows from the time-dependent Schrodinger equation with 
the Hamiltonian 



H = 



P 



1 muj 2 
2 



2m 2 4 
hfloa x cos(fci) cos(u;_l£) 



-2 1 

2gcos (tft)] x + -TiujQ(j z 



(1) 



Here the parameters a and q are proportional |?|] to the 
DC and AC voltages applied to the trap and a x ,cr z are 
the standard Pauli matrices. Moreover, we denote the 
frequencies of the AC field, the atomic transition, the 
Rabi frequency by w,cjo, and Qq, respectively. 

We introduce the dimensionless position x = kx, time 
t = Lui/2 and momentum p = -^-p. When we trans- 
form into an interaction picture with the unitary trans- 
formation exp(iu>LO' z t/2) the dimensionless Hamiltonian 
in rotating wave approximation reads 
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H 



4k ~ 1 1 

= -p 2 + - [a+ 2qcos(2t)} x 2 - kAa z 

mu! z 2 2 



kflo& x cos x. 
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Here we have defined the dimcnsionless Rabi frequency 
ft = ^o/ w an( i detuning A = (w^ — uq)/oj. 

The dynamics of the ion follows from the time- 
dependent Schrodinger equation 



««»=2r|*(t)) 



at 



for the state vector 



*(t)> = / dx[V g (x,t)\g) + V e (x,t)\e)]\x). 



(3) 



(4) 



order to investigate localization we start with this wave 
packet at the origin where the classical phase space is 
a stochastic sea @Jjl|. We calculate the time evolution 
using the split-operator method |2(J with a grid of 8192 
points. We control numerical errors using an adaptive 
time step-size algorithm 21 1. 

This numerical integration allows us to find the time 
dependence of the position probability distributions 
P g (x,t), P e (x,t) and P(x,t). Moreover, we calculate the 
corresponding momentum distributions where the wave 
functions $ g and \I/ e in position space are replaced by 
the wave functions in momentum space. 

We compare and contrast these quantum results to the 
dynamic s |B2| resulting from the classical equations of 
motion [|19|j23f| 



Here the effective Planck constant k = 2k 2 Ti/(muj) 
is consistent with the commutation relation [x,p] — 
2i 2 /(mw)[i,p] = 2k 2 /(muj)ih = ik. 

The state vector \^>) provides the probabilities 



or 



Pg{x,t)= \^g(x,t)\' 



P e { X ,t) = \^ e (x,t)\< 



(5) 



(6) 



to find the ion at the time t at position x given it is in the 
internal state \g) or |e), respectively. When we are not 
interested in the internal states the position probability 
distribution P(x, t) reads 



P(x,t) = \* g (x,t)\ 2 + \V e (x,t)\< 



(7) 



We simulate the effect of spontaneous emission by the 
quantum Monte-Carlo method 1 16 1^] using the effective 
non-hcrmitian Hamiltonian 



H, 



eff : 



H 



7 - - 
ik— o+o- 



x = p 

p = —(a + 2qcos(2t))x + fcf2 sin(x)ri (8) 

for the center-of-mass motion. These equations are 
driven by the Bloch equations p4| 



h = -2Ar 2 

r-2 = 2Ari + 2^o cos(a;)r3 
r3 = — 2Qq cos(a;)r2 



(9) 



describing the internal dynamics. Here r±, 7*2, and r$ de- 
note the in- and out-of-phase quadratures of the dipole 
moment, and the atomic inversion, respectively. 

For the comparison with the quantum mechanical re- 
sults we calculate 4096 trajectories starting from a clas- 
sical Gaussian ensemble centered initially at the origin 
and having the same widths in position and momentum 
as the quantum wave packet. In this way we obtain the 
classical phase space distribution from which we find by 
integration over x or p the classical momentum or posi- 
tion distributions P c i(p) or P c i(x). 



Here we introduced the atomic operators a+ = \e){g\ and 
<7_ = |<?)(e|, and 7 is the spontaneous decay rate scaled 
by U3. 

The moments of time when a spontaneous emission 
event takes place are chosen at random. Then the wave 
function is projected onto the ground state and renormal- 
ized, that is I*) — > The recoil p € [-k,k] 

is chosen randomly according to the probability distribu- 
tion fTJ 



of dipole radiation. When the results of single runs are 
averaged, we obtain the same result as predicted by a 
master equation |l(]{l7j. 

As an initial condition for the internal states we use 
the superposition (\g) + |e))/v2 and the the center-of- 
mass wave function is a Gaussian of width Aa; = \fk. In 



III. INFLUENCE OF THE DETUNING 

In the present section we study the influence of the 
detuning of the laser field with respect to the atomic 
transition on the dynamics of the system. Here we con- 
sider the following quantities of interest: the classical and 
quantum distributions, and their width as a function of 
time and detuning. 

In Fig. 1 we show the influence of the internal structure 
of the ion on the classical and on the quantum diffusion. 
Here we analyze for various detunings A the time depen- 
dence of the widths Aa; (left column) of the classical (thin 
lines) and quantum mechanical (thick lines) position dis- 
tributions P c i(x,t) and P[x,t) (right column). Whereas 
the left column emphasizes the time dependence of the 
widths, the right column shows position distributions av- 
eraged over the time interval [4507T, 5007r]. In order to 
bring out most clearly the influence of the detuning we 
have neglected spontaneous emission in this figure. 
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For a very small detuning the quantum width Ax lies 
well below the classical one. The latter increases as a 
function of time. In contrast, the quantum result dis- 
plays oscillations around a steady-state value. This sup- 
pression of the classical diffusion is due to dynamical lo- 
calization as discussed in Ref. ||. 

However, for slightly larger detunings the quantum 
curve is still oscillating but reaches partially above the 
classical one. For even larger detunings the quantum 
curve falls again below the classical one. Eventually for 
larger detunings it again goes partially above. Hence, for 
a fixed detuning, there exist time regimes in which the 
quantum diffusion is stronger than the classical one. 

This detuning dependence of the quantum diffusion 
manifests itself in the shape of the quantum mechani- 
cal position distribution. The latter consists of a sharp 
peak and a broad background. The detuning essentially 
controls the shape of the background: Depending on A 
we find distributions either of negative or positive cur- 
vature. In contrast, the detuning hardly influences the 
classical distributions shown on the right column by thin 
curves. We emphasize that the momentum distributions 
not shown here display almost identical behavior. 

To discuss the influence of the detuning on the local- 
ization length we show in Fig. 2 the widths Ax of the 
classical and the quantum distributions averaged from 
t = 2007T to 5007T. We note characteristic oscillations in 
the quantum mechanical curve as a function of the de- 
tuning A. These oscillations do not appear in the classi- 
cal curve which decreases monotonously for increasing A. 
This feature becomes clear when we recall that the width 
of the classical distribution at a given time is determined 
by the classical diffusion constant. The latter is propor- 
tional to the perturbing potential caused by the laser 
field. In the limitingcase of large detuning, the effective 
potential becomes |25[ proportional to the effective cou- 
pling constant SIq/ A. As A increases, the perturbation 
decreases and hence the diffusion is slower. 

An interesting domain is a small region around A = 0. 
Here the classical diffusion rate sharply decreases. The 
explanation of this feature follows from the set of Eqs. 
(^) and (||): The effective potential provided by the laser 
field in Eq. (||) is proportional to the factor r\{t). This 
quantity can take the maximum value of unity. For the 
resonant case, A = 0, we have r\ = and hence r\ is 
independent of time and determined by the initial condi- 
tion. The Bloch vector coordinates corresponding to the 
initial state (\g) + \e))/V2 are n(0) = l,r 2 (0) = and 
r 3(0) = 0. Hence r% = 1 results in the maximum driving 
strength of the center-of-mass motion. 

The curve for the quantum widths shows a distinctly 
different behaviour with a lot of structure. In particular, 
for A < 1 we note a well-developed resonance struc- 
ture and two sharp minima and maxima. The mini- 
mum around A = is related to the separable quan- 
tum dynamics in the two diabatic potentials = 
h\a + 2q cos(2£)] x 2 ± htto cos x. 



\g), |e) basis states to their superpositions |+) = (\g) + 
\e))/V2, and |-) = (\g) - |e))/\/2. In this basis the 
Hamiltonian Eq. (||) reads 



1 



1 



H = + - [a + 2qcos(2t)] x 2 
+fcSlo<T z cosi + kAa x , 



(10) 



and the coupling is now proportional to the detuning A. 
Hence we indeed find the potentials 

When we now increase the detuning, transitions be- 
tween the states |+) and |— ) contribute to the quantum 
diffusion. Therefore, the quantum diffusion rate quickly 
rises to the maximal classical one; this corresponds to the 
formation of the first maximum. 

The next minimum fits the condition A ps ui s , where 
u> s = 0.29 is the secular frequency of the trap. The ori- 
gin of this resonance becomes clear, when we write the 
Hamiltonian, Eq. (|^), in the interaction picture where the 
dynamics of the trap as well as the internal state energy 
are transformed away by unitary transformations ]2l| . In 
this picture the Hamiltonian reads 



oo oo oo 



#int(*) = E E E kJr +2k) ^ i - k ^ +A)t x 

n=0 k=-[n/2] l = ~oo 

a+\n)(n + 2k\ + h.c. (11) 

1 /2 

with the Lamb-Dicke parameter 77 = [k/{2uj r )] and 



(n,n-\-2k) <-> 



[zr?] - x 



(n + 2fc)! L 11 7r 



ir/2 



dt^^t)} 2 ^-^ 2 ^ 2 ^^ 2 ^)^- 2111 . (12) 

tt/2 

Here \n) denotes the n-th energy eigenstate of the time 
independent reference oscillator p7{ with frequency lo t . 
Moreover, we have made use of the Floquet solution 
e{t) = exp(io; s i)0(f), where 4>{t) = 4>(t + tt) is a peri- 
odic function. The solution e(t) obeys ^7]] the differential 
equation 



e+[a + 2qcos(2t)]e = 



(13) 



with e(0) = 1 and e(0) = iio r . 

From the term e Mi-ku s +A)t j n ^ e Hamiltonian 0) 
we expect resonance effects, when A satisfies the condi- 
tion I — ku) s + A = 0. For I = 0, k = 1 we have A = u> s , 
which corresponds to two-phonon transitions. This reso- 
nance suggests an enhancement of the quantum diffusion 
rate on resonance. Instead, we have a deep minimum. 

In order to understand this counter-intuitive behav- 
ior we calculate the characteristic frequencies w | n, ™ +2 ' c ) 
for I — and k = 1, 2 at different vibrational quantum 

numbers n. From the inset of Fig. 2 we recognize that 

I, .,(«,n+2). 



has a deep minimum around n — 10, which 
To bring this out most clearly we transform from the explains suppression of the diffusion over the vibrational 
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states. As shown by the curve in the inset for \cOq n | 
there is no minimum just a decay. This causes the quan- 
tum diffusion to take on the classical value, which ex- 
plains the formation of the second maximum [ p8[ . 

We conclude this section by noting that the phe- 
nomenon of oscillations in the width of distributions also 
appears in another quantum system showing classical 
chaos and quantum localization: an atom moving in a 
phase-modulated standing wave [p|JlO|. However in this 
case the oscillations are of classical origin and appear 
both in the classical and quantum widths. 

IV. INFLUENCE OF SPONTANEOUS EMISSION 

We now consider the effect of spontaneous emission 
using the quantum Monte-Carlo technique (l^JT^]. The 
purpose of our investigations is twofold. First, we want 
to show that indeed the phenomena discussed so far are 
quantum interference effects. They are therefore sensitive 
to decoherence such as spontaneous emission. The latter 
predominantely occures when the laser is tuned close to 
resonance with the atomic transition. Second, we want 
to show that in the far-detuned limit the effect of deco- 
herence is negligible and the phenomenon of dynamical 
localization survives. 

In Fig. 3 we show the results of our simulations with 
a realistic rate of spontaneous emission corresponding to 
the decay rate of the S — > P transition at 19.4 MHz 
m 9 Be+ ]2|. We compare them to the classical and to 
the quantum result without noise. The loss of the co- 
herence causes destruction of localization. In this case 
the widths of the quantum mechanical position and mo- 
mentum distributions are larger than the classical ones. 
This additional diffusion is caused by the random recoil 
kicks following each spontaneous emission event. When 
we neglect the recoil — which of course is not realistic — 
the quantum curve follows the classical one. The position 
and momentum distributions are of classical type, that 
is on the logarithmic scale they are polynomial curves. 
Note that the pattern of the standing wave appears on 
top of the position distributions. 

In Fig. 4 we show the results for A = 1000, correspond- 
ing to a detuning of 10 GHz, a value which was mentioned 
in Ref. [^| . The small oscillations in the quantum widths 
Aa; and Ap are destroyed, but the main phenomenon, the 
substantial quantum suppression of classical diffusion is 
still visible. 

It is interesting to note that in a single realization of 
the dynamics the coherence of the two-level superpo- 
sition is completely destroyed by a single spontaneous 
emission. However, it affects only slightly the motional 
coherence because the population of the excited state is 
very low. Furthermore, the small difference between the 
results with and without spontaneous emission in Fig. 4 
is of the same order of magnitude that was found for the 
problem of an atom in a phase modulated standing wave 




In the position distribution of the ground state, there 
is more probability in the classical-like background than 
in the case of no spontaneous emission. We have found 
that this background is very slowly growing, and no real 
steady state can be reached. However, if we increase the 
detuning, the real steady state is approached asymptot- 
ically. Since dynamical localization appears for a large 
range of parameters || , there is a lot of room for optimiz- 
ing the parameters. Thus the phenomenon of dynamical 
localization in a Paul trap could be observed experimen- 
tally. 

In order to observe not only dynamical localization but 
also the oscillations in the localization length discussed in 
the previous section, one should consider a configuration 
where decoherence is weaker even for small detunings, 
for example a different set of parameters in the present 
system, or a different system such as a Raman transition 
between two ground states. 

V. CONCLUSIONS 

There are basically two ways of controlling dynamical 
localization: (i) through the classical diffusion and (ii) 
through the quantized nature of the variable showing lo- 
calization. This stands out most clearly in the estimate 
/ r~> D/k 2 for the localization length. Here D is the clas- 
sical diffusion coefficient determined by the perturbing 
potential. The coupling to this potential is called the 
control parameter, since it is a direct way to control the 
localization length. The scaled Planck constant k de- 
scribes how important the quantization of the system is 
with respect to the perturbation. 

In this paper we have shown that, when the perturbing 
potential has a quantum character as well, the relation 
I ~ D/k 2 is not exactly true anymore; we observe os- 
cillations in the localization length which do not appear 
in the classical diffusion rate. The parameter determin- 
ing the oscillations could be called the quantum control 
parameter. 
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FIG. 1. Classical and quantum dynamics of the cen- 
ter-of-mass motion of a two-level ion stored in a Paul trap 
and interacting with a classical standing light wave. We 
show the influence of the detuning A between the frequen- 
cies of the atomic transition and the light on the dynam- 
ics using probability distributions P(x) in position averaged 
over the time period [450-7T, 50O7r] (right column), and widths 
Aa; = [(a; 2 ) — (a;) 2 ] 1 ^ 2 (left column) as a function of time. 
For a very small detuning (top) the quantum mechanical 
width (thick line) lies below the classical width (thin line). 
Whereas the width of the classical position distribution in- 
creases monotonously as a function of time the correspond- 
ing quantum result displays oscillations around a steady-state 
value. This suppression of the diffusion is due to dynamical 
localization as discussed in Ref . g . For slightly larger detun- 
ings the quantum curve is still oscillating but goes above the 
classical one. For even larger detunings the quantum curve 
falls again below the classical one and the diffusion is sup- 
pressed by dynamical localization. However, for larger detun- 
ings (bottom) the quantum curve again starts to move across 
the classical one. We note that only the quantum curves 
show a strong dependence on the detuning. This dependence 
also manifests itself in the position distributions. Indeed A 
influences the shapes of the quantum, but not of the clas- 
sical distributions. Here and in Fig. 2 we do not account 
for spontaneous emission. We have used the scaled Rabi fre- 
quency flo — 2.24 of the standing wave and the trap param- 
eters a = 0, q = 0.4 corresponding to the secular frequency 
ui s — 0.29. The effective Planck constant is k — 0.29. 

FIG. 2. Dependence of the classical (dashed line) and quan- 
tum mechanical (solid line) variance Aa; in position on the de- 
tuning A. For each value of A we have averaged the variances 
over the time window [2007T, 5007r]. Whereas the classical 
curve decays monotonously with a characteristic maximum 
on resonance, the quantum curve shows striking oscillations. 
The inset shows the amplitudes |tjQ n,n+2 '| and |tjg n ' n+4 ' de- 
termining the strengths of the 2-phonon and 4-phonon tran- 
sitions in their dependences on the vibrational quantum num- 
ber n. Whereas |a;Q n ' n+2 '| has a minimum, \uj^ l ' n+A ^\ decays 
monotonously. This leads to a minimum or a maximum in 
Aa; for A = ui s or A = 2uj s , respectively. The squares mark 
the values of Aa; at the detunings used in Fig. 1. Here and in 
Fig. 3 the parameters are as in Fig. 1. 



FIG. 3. Influence of spontaneous emission on dynamical lo- 
calization for A = 0, that is on atomic resonance: We show 
the time dependence of the widths Aa; and Ap in position and 
momentum (left column) and the probability distributions 
P(x) and P(p) averaged over the time interval [2007T, 2507r] 
(right column). Fat, thin and thick curves indicate these 
quantities in the presence and absence of spontaneous emis- 
sion and for the classical case, respectively. We note that 
spontaneous emission destroys localization and due to the re- 
coil the quantum widths become even larger than the clas- 
sical ones. In the absence of spontaneous emission the posi- 
tion distributions of the ground and excited state (jagged thin 
curve) are sharply peaked. Moreover they are identical due 
to the fact, that the initial internal state is the superposition 
(| e) + \g) ) /v2- Spontaneous emission destroys the localization 
peak and enhances the population in the ground (upper fat 
curve) at the expense of the ions in the excited state (lower fat 
curve). In this case the quantum distributions in position and 
momentum are even broader than the corresponding classical 
ones (thick lines). The pattern of the standing wave reflects 
itself in the position distributions (upper right corner). It is 
more pronounced in the quantum mechanical curves than in 
the classical one. The results in the case of spontaneous emis- 
sion have been obtained by averaging over a sufficient number 
of runs, in this case 79 runs. Here and in Fig. 4 the scaled 
rate of spontaneous emission is 7 = 2. 

FIG. 4. Influence of spontaneous emission on dynamical lo- 
calization for A = 1000, that is far off-resonance. As in Fig. 
3 we show the widths in position and momentum as a func- 
tion of time and the corresponding probability distributions 
averaged over the time interval [2007T, 2507r]. Since we start 
from the ground state of the ion the strong detuning prevents 
the population of the excited state. We therefore only display 
the ground state population (fat line). We emphasize that 
this curve differs only slightly from the corresponding curve 
in the absence of spontaneous emission (thin line). Sponta- 
neous emission does not destroy the localization. Due to the 
small amount of spontaneous emission an average over 49 runs 
was sufficient. In order to make contact with Ref. we used 
the parameters Oo = 94.69 and k — 0.0725. 
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